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Abstract. We improve the recent result of 2 proving a one-sided threshold condi- 
tion which leads to finite-time breakdown of the Euler-Poisson equations in arbitrary 
dimension n. 



1. Introduction 
The pressure-less Euler-Poisson (EP) equations in dimension n > 1 are 

(1.1a) pt + div (pu) = 

(1.1b) ut + u- Vu = kVA'^{p- c), 

governing the unknown density p = p{t,x) : IR+ x i— ]R_^ and velocity u = u(t,x) : 
]R+ X M" I— > M" subject to initial conditions p(0,x) = po{x) and u(0,x) = Uo(x). They 
involve two constants: a fixed background state c > such that j {p — c)dx = 0, and the 
constant k which parameterizes the repulsive k > or attractive k < forcing, caused 
by the Poisson potential A~^(p — c). The EP system appears in numerous applications 
including semiconductor, plasma physics {k > 0) and collapse of interstellar cloud [k < 0). 
This paper is restricted to the attractive case, k < 0. For simplicity, we set c = 1, 



k = —1 in (1.1a), (1.1b) to arrive at the unit-free EP system, 

(1.2a) Pt + div (pu) = 0, 

(1.2b) ut + u- Vu = -VA"^(p- 1). 

All discussion in this paper remains valid for EP system with physical parameters c > 0, 



A; < upon a simple rescaling argument — see Corollary 1.1 below. 

We are concerned here with the persistence of regularity for solutions of the attrac- 
tive EP system. Our main theorem reveals a pointwise criterion on the initial data, a 
so-called critical threshold criterion [5l El HO], that leads to finite time blow-up of Vu. 
It is a sharp, nonlinear quantification of balance between divn and p, two competing 
mechanisms that dictate the regularity of EP fiows. Our result also stands out as a 
generalization of several existing results O |2l [9l [10] for which further discussion is given 
after the Main Theorem and its corollary. 
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Main Theorem 1.1. Consider the n-dimensional, attractive Euler-Poisson system (1.2a), 



(1.2b) subject to initial data po, uq. Then, the solution will lose regularity at a finite 
time t = tc < oo, if there exists a non-vacuum initial state po{x) > with vanishing initial 
vorticity, V x uq{x) = 0, such that the following sup-critical condition is fulfilled. 



(1.3a) 

where 



divuQ{x) < sgn{pQ{x) — l)^/nF{pQ{x)), 



(1.3b) F{p) 
In particular, min^ div u(t, x) 



n 



np- 



n 



2/n 



1 -p + plnp, 
—oo and max^. p(t, x) 



n^2, 
n = 2. 

■> cxD as t ] tc- 



Proof. Combine Lemma 3.1 and Lemma 4.2 below, while noting that the curve 

div u = sgn{p — 1) a/ nF{p), 
is the separatrix along the boundary of the blow-up region f2 = fii U defined in (4.3) 



and illustrated in Figure 4.1 



□ 



We note by passing that, by classical arguments, the force-free Euler system Ut-fu- Vu = 
exhibits finite time blow-up if and only if there exists at least one negative eigenvalue of 
Vuo(x). In the above theorem, however, finite-time blow-up can occur solely depending 
on the initial profile of div uq and po regardless of individual eigenvalues of Vuq. 

We also note that, by rescaling p to p/c, x to \/—kcx and t to \/—kct, the Main 



Theorem immediately applies to the EP system (1.1a), (1.1b) with physical parameters. 
Since the EP system with k < models the collapse of interstellar cloud, the following 
corollary reveals a pointwise condition for mass concentration, p oo, which interestingly 
preludes the birth of new stars. 



Corollary 1.1. Consider the Euler-Poisson system (1.1a), (1.1b) with c > 0, A; < 
subject to initial data po, uq. Then, the solution will lose regularity at a finite time 
tc < oo, if there exists a non-vacuum initial state po{x) > with vanishing initial vorticity, 
V X uo(x) = 0, such that the sup-critical condition is fulfilled. 



;i-4) 



div uo{x) < sgn{po{x) — c)\ —nkcF 



where F{-) is given in (1.3b). In particular, min^. div u(t, x) 
as t t tc- 



-oo and max^ p{t, x) oo 



The concept of Critical Threshold and associated methodology is originated and de- 
veloped in a series of papers by Engelberg, Liu and Tadmor [3] , Liu and Tadmor [TUl |H] 
and more. It first appears in [3] regarding pointwise criteria for solution regularity 
of ID EP system. The key argument in that paper is based on the convective derivative 
along particle paths ' = dt + u ■ V. It makes possible to obtain a 2-by-2 ODE system 
for Ux and p along particle paths — the so-called Lagrangian formulation. Phase plane 
analysis is then employed to study the finiteness of the ODE solutions and therefore 
regularity of the PDE solution. Similar results stay valid for Euler-Poisson systems with 
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geometric symmetry in higher dimensions (HE]. To treat genuinely muhi-D cases, Liu 
and Tadmor introduce in [8] the method of spectral dynamics which relies on the ODE 
system governing eigenvalues of 

M := Vu, 

which is the velocity gradient matrix, along particle paths. They identify if-and-only- 
if, pointwise conditions for global existence of solutions to restricted Euler-Poisson 
systems. Chae and Tadmor [2] further extend the Critical Threshold argument to multi- 



D full Euler-Poisson systems (1.2a), (1.2b) with attractive forcing A; < 0. Their result. 



however, offers a blow -up r egion V x uq = 0, div uq < —^—nkc which is only a subset of 



the blow-up region in (1.4). This subset is to the left of the solid line d < d^ := —\J—nkc 
depicted in figure 4.1 Finally, a recent paper by Tadmor and Wei [22] reveals the critical 
threshold phenomena in ID Euler-Poisson system with pressure. 

When tracking other results on well-posedness of Euler-Poisson equations, we find them 
commonly relying on (the vast family of) energy methods and thus fundamentally differ 
from our pointwise results obtained via the Lagrangian approach. With repulsive force 

> 0, we refer to [?', 'S] for global existence of classical solutions with small data and 
[19] for nonexistence of global solutions. With attractive force /c < 0, see [I5j for local 
regularity of classical solutions and [161 E] foi' nonexistence results. Discussion on weak 
solutions of Euler-Poisson systems can be found in e.g. [251 HH ED]. We also refer to 
[HI m [121 [in 121] and references therein for steady-state solutions. Study of Euler-Poisson 
system with damping relaxation can be found in e.g. [23l |2H [13] . 

The rest of this paper is organized as following. In Section |2| we follow the idea of 
[2] to derive along particle paths an ODE system governing the dynamics of eigenvalues 

for S := -(M + M^). This is a variation of the spectral dynamics for M introduced in 



[H]. We then derive in Section [3] a closed 2x2 ODE system (3.1) at the cost of turning 
one equation into inequality. By the comparison principle, this inequality is in favor of 
blow-up. Thus, with the inequality sign being replaced with equality sign, a modified 
ODE system is used to yield sub-solutions and to study blow-up scenario for the original 
system. Section |4| devoted to the modified system, reveals the Critical Threshold for 
such a system. Consequently, a pointwise blow-up condition for the original system is 
identified. 



2. Spectral dynamics 

We examine the gradient matrix M = Vu and its symmetric part, S = ^ (Vu + (Vu)"'") . 
Both matrices are used to study the spectral dynamics of Euler systems (see e.g. [8] for M 
and [2] for S). The relation between the spectra of M and S is described in the following. 

Proposition 2.1. Let {A^} denote the eigenvalues of M and {Xs} for S. Then 

(2.1) Aj\/ = As = div u. 

Am 

(2.2) E^M = E4-^M^- 

Am Ag 
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Here, u is the vorticity vector which consists of the off- diagonal entries of A :- 

i(vu-(vun. 

Proof. Use identity M = S + A and the skew-symmetry of A, 

Am = tr{M) = tr{S + A) = tr{S) = ^ A5. 

Squaring the last identity we have = S'^ + A"^ + AS + SA and therefore, 
^ XIj = tr{M^) = tr{S^ + A^ + AS + SA) = J2 + tr{A^)- 



Note that AS + SA is skew-symmetric and thus traceless. A simple calculation yields 
tr{A^) = □ 

Following [8j, we turn to study the dynamics of M along particle paths. Take the 
gradient of (1.2b) to find 

(2.3) M' + M^ = Mt + u ■ VM + = -R{p - 1), 

where R stands for the Riesz matrix, R = {Rij} ■= {dxiXj^"^}- 

The trace of (2.3) then yields that the divergence, d := tr{M), is governed by 

^' = -E^M-(p-i), 

Am 



and in view of (2.2) 
(2.4) 



d' 



UJ\ 



(p-1)- 



We now make the first observation regarding the invariance of the vorticity 00: taking 



the skew-symmetric part of the M-equation (2.3), 
(2.5) A' + AS + SA = 0. 



It follows that if the initial vorticity vanishes, uj{x) V x u(x) = 0, then by (2.5), 
a; H- V X u vanishes along the particle path which emanates from x. This allows us to 



decouple the vorticity and divergence dynamics, and (2.4) implies 



(2.6) 



d' 



J]A|-(p-l), Vxu = 0. 



As 



Finally, we use Cauchy-Schwartz A| < — 



2 X 

A5 ) = —d"^ and the fact that all Xs 
/ n 

are real (due to the symmetry of S), to deduce the inequality, 



(2.7a) 



d' < d^-(p- 1) 

n 



This, together with the mass equation (1.2a) which can be written along particle path 
(2.7b) p = -dp, 

give us the desired closed system which governs (p, d) along particle paths. 
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Remark 2.1. The approach pursued in this paper will he based on the inequality (2.7a) 
and is therefore limited to derivation of a finite time breakdown. To argue the global 



trace ^ A|. In the two-dimensional case, for example, one can use ^ A| = d'^ /2 + rf /2 



regularity, one needs to study the underlying equality {2.0^ , and to this end, to study the 
trace ^ A|. In the tw 
to replace (2.7a) with 

d' 



1 



1 



7] :-- 



As,2 — A 



S,l- 



In this framework, global 2D regularity is dictated by the dynamics of the spectral gap, 
V = -^3,2 ~ ^s,i! which in turn requires the dynamics of the Riesz transform R{p — 1). 

3. A COMPARISON PRINCIPLE WITH A MAJORANT SYSTEM 



The blow-up analysis, driven by the inequalities (2.7), 

(3.1a) d' < --d^ - (p-l), 

n 

(3.1b) p' = -dp. 

is carried out by standard comparison with the majorant system 

(3.2a) e' = --e2-(C-l), 

n 

(3.2b) C' = -e(. 

The following proposition guarantees the monotonicity of the solution operator associ- 
ated with (3.1 ). 



Lemma 3.1. The following monotone relation between system (3.1) and system (3.2) is 
invariant forward in time. 



(3.3) 



rf(0) < e(0) 

< C(0) < p(0) 



implies 



d{t) < e{t) 
< at) < pit) 



for t > 0, 



as long as all solutions remain finite on time interval [0,t]. 



Proof Invariance of positivity of C is a direct consequence of (3.2b) and finiteness of e. 
The rest can be proved by contradiction. Suppose ti is the earliest time when (3.3) is 
violated. Then, 

(3.4) C{ti) = C(0) exp e{t)dt^ < p(0) exp d{t)dt^ = p(ti). 

Therefore, we are left with only one possibility e{ti) = d{ti). However, subtracting (3.1a) 

1 



from (3.2a), 
(3.5) 



{e-dy>—ie'-d')-{C-p). 



n 



Setting t = ti in the above inequality, we find that 

LHS of ^ = (e(ti) - d{ti)y < 0, 

since e(t) — d(t) > for all t < ti; but this contradicts the positivity of the expression on 
the right of (3.5), for by (3.4) 

RHS of K§ = - [C{ti) - p{ti)] > 0. 
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□ 



In the next section, we employ phase plane analysis on the modified system (3.2). When 



translated in terms of the original system (3.1 ), however, such analysis can only yield blow- 



up results and is insufficient for global existence results. In other words, estimate (3.3) is 
only useful for proving d \ — oo, the key mechanism for blow-up of solutions. 

4. Stability analysis of the majorant system 



We shall prove the blow-up of the majorant system (3.2), e{t) — oo as t | tc, which 
in turn, by lemma 3.1 implies d{t) — oo. Abusing notations, we express the majorant 
system in terms of the original variables (e, Q ^— (rf, p): 

(4.1a) d' = --d'^ - (p 



(4.1b) 



d' 
I 

P 



-d^ 
n 

-dp. 



The (in-) stability analysis of (4.1) hinges on the path invariants of this system. To 
this end, we use the same g-transformation employed in [HI (TU]: setting q := d^ and 
differentiate along the path {{t, X{a,t)) \ Xt{a,t) = u{t, X{a,t)), X{a,0) = a}, we find 



dq 
dp 



2d- 



1-- — 
p' np 



g + 2 1 



which yields 



d 

dp 



qp 



=2(l-p-V 



Upon integration, we arrive at the following key observation. 



Lemma 4.1. The majorant system (4-1) is equipped with the path invariant, 



/(rf(t),p(t)) = /(c/o,Po), 
along each path {t,x{t)) initiated with a non-vacuum state {do, po > 0). Here, 



(4.2) 



I{d, p) := d'^p 



(1 — r ^)r " dr = p " (d^ — nF{p)^ 



where F{-) is specified in (1.3b). 



It is simple calculation to show that the majorant system (4.1) admits three distinct 



critical points (see figure 4.1) 



(d*,p*):=(0,l), (d±,p±):=(±y^,0). 

and that (0, 1) is a saddle point, {—y/n, 0) a nodal source and {y/n, 0) a nodal sink. The 
separatrix is given by the zero level set I{d,p) = 0. Moreover, the right branch of the 
separatrix, d = ^JnF{p) connects critical points (0, 1) and {\/n, 0) while the left branch, 
d = —\JnF{p) connects (0, 1) and {—^/n, 0). 



By inspection of the phase plane in figure 4A, we postulate the following invariant 
region of finite-time blow-up for the modified system ( |4.1 ), 

(4.3a) 



VL = VLiUVL2 = {{d,p) I d < sgn{p - l)^JnF{p)] 
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where 

(4.3b) 
(4.3c) 



Qi := {{d,p) I I{d,p) > and d <0 and p > 0}, 
^2 := {{d, p) \ I{d, p) < and p > 1}. 



t 
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Figure 4.1. Phase plane of (4.1) with blow-up region U which ex- 



tends the Chae-Tadmor region \1\ d < d 



Lemma 4.2. Consider the modified system (4-1), equipped with initial data {do,po). If 



ido,Po) £ then divu — > —oo and p ^ oo at a finite time. 



Proof. We begin by recalling (1.3b), consult (4.2) 



F{p) 



n 



-pn 



(1 — r ^)r " dr. 



Clearly, -F(l) = F'(l) = and a simple calculation shows that F"ip) = —p^ ^, which 

n 

implies that F{p) is a strictly convex function of positive p and attains its only minimum 

at p = 1, 



(4.4) 



F{p) > F(l) = 0. 



We shall also utilize the invariance of (4.2) 

(4.5) d^ - nF{p) = pn/o, /q = /(^o, po)- 

We now turn to discuss the two possible blow-up scenarios, depending whether the 



initial data ((io,Po) belong to the blow-up regions fii or given in (4.3). 
Case #1. Assume that {do, po) G so that the invariant / remains a positive constant 

/ > 0. 



8 



BIN CHENG AND EITAN TADMOR 



In this case, d remains negative 
Fip)- 



'or otherwise, setting d = in (4.5) would resuh in 



pnl /n < 0, violating (4.4). Thus, (4.5) and (4.4) yield an upper bound, 

d < —p"^/l. 



Then, by (4.1b), we have a Riccati type of equation p' > -\/7p^^" for which the solution 
exhibits blow-up p — *• +oo and the divergence d = divu approaches —oo at a finite time 



due to (4.5) 



Case #2. Assume that {do, po) G ^2 so that the invariant / remains a negative constant 

/ < 0. 

1 remains positive, for otherwise setting p 



In this case, p 

F{1) = [d"^ — I)/n > in contradiction to (4.4). Now, for p > 1 we have 
F{P) = -p''- r(l--] 4u;dr < ^-p^'-{p - 1 



1 in (4.5) would result in 



n 



n 



-P 



This together with (4.5) yield 

%2/"(p-l)>F(p) 

- 1 > — //2 follows. Thus, by (4.1a), we end up with a Riccati 

d^ 



n 



(d'-p^/^l) > --p^/^I 



and the lower bound, p 
type of equation 



d' < 



I 

^ ^ 2' 

Since the invariant / remains a negative constant, the solution exhibits blow-up d = 
div u ^ — oo at a finite time even if initially d^ > 0. The density p also approaches oo in 
finite time due to (4.5). □ 



The last step of proving the Main Theorem is just to combine the comparison principle 



in Lemma 3.1 with the above lemma. We notice that fl is an open set and thus given any 



initial data (rfo, Po) £ ^ for the original system, we can always find e > and initial data 
{do + e, po — e) G Q for the modified system. This latter initial data will lead to finite time 
blow-up of the modified system and therefore, by lemma 3.1, initial data {do,po) G Q will 
lead to finite time blow-up of the original system. 
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